Lattice based equation of state and transverse momentum spectra of identified 
particles in ideal and viscous hydrodynamics 
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Assuming that in Au+Au collisions, a baryon free fluid is produced, transverse momentum spectra 
of identified particles (n, K, p and 0), in evolution of ideal and viscous fluid is studied. Hydro- 
dynamic evolution is governed by a lattice based equation of state (EOS), where the confinement- 
deconfinement transition is a cross-over at T co =196 MeV. Ideal or viscous fluid was initialised to 
reproduce <f> meson multiplicity in 0-5% Au+Au collisions. Ideal or minimally viscous (jy/s=0.08) 
fluid evolution reasonably well explain the transverse momentum spectra of pion's, kaon's, and <j> 
meson's in central and mid-central Au+Au collisions. Description to the data is much poorer in 
viscous fluid evolution with n/s >0.12. The model however under estimate proton production by a 
factor ~ 2. 

PACS numbers: 47.75+f, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 



Relativistic hydrodynamics provides a convenient tool 
to analyse relativistic heavy ion collision data. It is as- 
sumed that a fireball is created in the collisions. Con- 
stituents of the fireball collide frequently to establish lo- 
cal thermal equilibrium sufficiently fast and after a cer- 
tain time Tj, hydrodynamics become applicable. If the 
macroscopic properties of the fluid e.g. energy density, 
pressure, velocity etc. are known at the equilibration 
time Ti, the relativistic hydrodynamic equations can be 
solved to give the space-time evolution of the fireball till a 
given freeze-out condition such that interaction between 
the constituents is too weak to continue the evolution. 
Using suitable algorithm (e.g. Cooper-Frye) information 
at the freeze-out can be converted into particle spectra 
and can be directly compared with experimental data. 
Thus, hydrodynamics, in an indirect way, can character- 
ize the initial condition of the medium produced in heavy 
ion collisions. Hydrodynamics equations are closed only 
with an equation of state (EOS). It is one of the most 
important inputs of a hydrodynamic model. Through 
this input macroscopic hydrodynamic models make con- 
tact with the microscopic world and one can investigate 
the possibility of phase transition in the medium. Most 
of the hydrodynamical calculations are performed with 
EOS with a 1st order phase transition. For example, 
in [l|, 1st order EOS, EOS-Q was used to successfully 
analyse a host of experimental data in Au+Au collisions 
at RHIC. In EOS-Q, the deconfined or the Quark-Gluon 
Plasma (QGP) phase is modeled by a bag equation of 
state of non-interacting quarks and gluons, the confined 
or the hadronic phase is modeled by a non-interacting 
gas of hadronic resonances. Ideal hydrodynamics analy- 
sis with EOS-Q indicate that in central Au+Au collisions, 
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at the equilibration time Tj ~ 0.6 fm, central energy den- 
sity of the QGP fluid is e, «30 GeV/f m - 3 [1]. However, 
lattice simulations 0, Q indicate that the confinement- 
deconfinement transition is neither a 1st nor a 2nd order 
phase transition, rather a cross-over at T co —196 MeV. It 
is then important that lattice based EOS are used in hy- 
drodynamic analysis of RHIC data, more so when we are 
trying to verify the lattice prediction about confinement- 
deconfinement transition. 

One also note that ideal hydrodynamic predictions for 
the initial energy density or temperature of the fluid pro- 
duced in Au+Au collisions is not creditable as dissipative 
effects are not included. In hydrodynamics, initial energy 
density or temperature of the fluid is obtained by fitting 
experimental data on particle production, e.g. pion mul- 
tiplicity, pt spectra etc., which essentially measure the 
final state entropy. Unlike in ideal fluid evolution, where 
initial and final state entropy remains the same, in vis- 
cous fluid evolution entropy is generated. Consequently, 
to produce a fixed final state entropy, viscous fluid re- 
quire less initial entropy density (or energy density) than 
an ideal fluid. 

Recently, we have constructed an lattice based EOS 
and use it to explain the <fi meson production in Au+Au 
collisions at RHIC [||. Recent lattice simulation results 
were parameterised to obtain EOS of the deconfined 
state. The confined part of the EOS corresponds to that 
of a hadronic resonance gas with all the resonances with 
mass m res <2.5 GeV. The deconfined and the confined 
part of the EOS were smoothly joined at a cross-over 
temperature T co —196 MeV. In [J] it was shown that the 
lattice based EOS reasonably well explain the centrality 
dependence of <j> mesons multiplicity (dN^/dy), mean pt 
((pj)) an d integrated elliptic flow (v^)- From a simul- 
taneous fit to dN^/dy, (pj,) and v% an estimate of the 
shear viscosity to entropy ratio was obtained, f7/s=0.07 ± 
0.03 ± 0.14, the first uncertainty is due to uncertainty in 
STAR measurements, the 2nd one is due to uncertain ini- 
tial condition e.g. initial time varying between Ti=0.2-0.6 
fm), freeze-out temperature varying between 7V=130- 
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FIG. 1: Viscous fluid (r;/s=0.08) simulation for temporal 
evolution of momentum anisotropy in b~7 fm Au+Au col- 
lision at RHIC. The solid line is the simulation result from 
VISH2+1 [f| and the dashed line is the simulation result from 
AZHYDRO-KOLKATA. 



FIG. 2: Black circles are lattice simulation [2] for entropy 
density. The black line is the parametric representation to the 
combination of lattice and hadron resonance gas EOS joined 
smoothly at cross over temperature T CO =196.0 MeV. 



150 MeV, initial velocity (v r = tanh(ar), a=0.0-0.06), 
inaccuracy in hydrodynamic evolution code etc. 

Strange meson <fi constitute only a very small fraction 
of the total particles produced in Au+Au collisions. Par- 
ticle production is dominated by pions, kaons and pro- 
tons etc. For example, in central Au+Au collisions, pions 
constitute nearly ^80% of the total particle yield, Kaons 
~ 13% and protons ~ 5%. <f) mesons contribute ~ 2% to 
the total yield. It is then important to inquire whether 
or not hydrodynamics with the lattice based EOS is con- 
sistent with the experimental data on other particles, e.g. 
pion, kaon, proton etc. The estimate of viscosity as ob- 
tained in [4( will not be creditable unless the model also 
reproduces bulk of the particles, i.e. w, K, proton etc. In 
the present paper, with the same parameters as in [ij , we 
have analysed transverse momentum spectra of identified 
particles, e.g. n, K, proton and </> in Au+Au collisions 
over a wide range (0-60%) of collision centrality. In cen- 
tral and mid central collisions, hydrodynamic evolution 
of minimally viscous fluid best explain the data. Nearly 
equilvalent description is also obtained in ideal fluid evo- 
lution. However, description to the data in evolution 
with viscosity rj/s > 0.12 is considerably poor than that 
in ideal or minimally viscous fluid. 

The paper is organised as follows: in section [Hi we 
briefly describe the hydrodynamical equations used to 
compute the evolution of ideal and viscous fluid. Con- 
struction of the lattice based equation of state is dis- 
cussed in section III Bl Simulation results are discussed 
in section III CI Summary and conclusions are given in 
section QV] 



II. HYDRODYNAMICAL EQUATIONS, 
EQUATION OF STATE AND INITIAL 
CONDITIONS 

A. hydrodynamic equations 

In Israel- Stewart's theory of 2nd order dissipative hy- 
drodynamics, for a baryon free fluid, and neglecting bulk 
viscosity and heat conduction, the space-time evolution 
of a relativistic fluid is obtained by solving, 

= 0, (1) 

DTT^ = - — (7T^-27 ? V < ' 1 U I/> ) 

- [<it'V A + u u ti ijX ]Du x . (2) 

EqQ] is the conservation equation for the energy- 
momentum tensor, T^ v = (e + p)u^u v — pg^ v + 7r A ' !/ , 
e, p and u being the energy density, pressure and fluid 
velocity respectively, tt^ is the shear stress tensor. Eqj2] 
is the relaxation equation for the shear stress tensor . 
In Eql2l D — u^d^ is the convective time derivative, 
V<»u v> = ±(V»u v + V^ 1 ) - \{d.u){g^ v - u^u v ) is a 
symmetric traceless tensor. r\ is the shear viscosity and t„ 
is the relaxation time. It may be mentioned that in a con- 
formally symmetric fluid relaxation equation can contain 
additional terms [5]. Assuming boost-invariance, EqsQ] 
and [2] are solved in (r = \lt 2 — z 2 , x, y,r) s — | In co- 
ordinates, with a code '"AZHYDRO-KOLKATA" 1 , de- 
veloped at the Cyclotron Centre, Kolkata. Details of 
the code can be found in As shown in FigQ] for 

similar initial conditions, within 10%, the code repro- 
duces the temporal evolution of momentum anisotropy 
£ p = <r^+r««> °f a QGP nm d as calculated by Song 
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FIG. 3: Ideal hydrodynamic predictions for pr spectra 
of - in 0-5% Au+Au collisions are compared with 

PHENIX data [Tq| . Thermal and decay pions are shown sep- 
arately. 

and Heinz @. 

B. Equation of state 

Equation of state (EOS) is one of the most impor- 
tant inputs of a hydrodynamic model. Through this 
input macroscopic hydrodynamic models make contact 
with the microscopic world. Most of the hydrodynami- 
cal calculations are performed with EOS with a 1st order 
phase transition. Huovinen [l8| reported an 'ideal' hy- 
drodynamic simulation with a cross-over transition. He 
concluded that the experimental data (e.g. elliptic flow 
of proton or antiproton) are better explained with EOS 
with 1st order phase transition than with EOS with 2nd 
order phase transition. Huovinen [lj| used the 'thermal 
quasiparticle model' [lj| to obtain EOS for the decon- 
fined phase. For the confined phase he used the hadronic 
resonance gas model. 

Recently, Cheng et al [2| presented high statistics lat- 
tice QCD results for the bulk thermodynamic observ- 
ables, e.g. pressure, energy density, entropy density etc. 
The simulations were performed with two light quarks 
and a heavy strange quark. The quarks masses are 'al- 
most' physical, and corresponding pion mass is m n ~ 
220 MeV. The strange quark mass was adjusted to phys- 
ical value rriK ~503 MeV. In FigJ5J we have shown the 
simulation result for the entropy density Q. We have 
parameterise the entropy density as, 

7 L = a+[0 + 1 T][l + tanh^^}, (3) 
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FIG. 4: (color online) PHENIX data [ijj on the central- 
ity dependence of 7r+ + 7r in 0-60% Au+Au collisions. For 
better visibility experimental data beyond 0-5% are divided 
by a factor of 10 respectively. The solid, dashed, medium 
dashed and long dashed lines are hydrodynamic predictions 
from ideal(?7/s=0.0) and viscous (??/s=0.08,0.12,0.16) fluid re- 
spectively. 



In Fig. 2, the solid curve is a parameterisation with 
q=0.64, ^=6.93,7=0.55, T c =196MeV,AT=0.1T c . From 
the parametric form of the entropy density, pressure and 
energy density can be obtained using the thermodynamic 
relations, 



p(T) = [ s{T>)dT> (4) 
Jo 

e{T) = Ts-p. (5) 

We complement the lattice EOS 01 by a hadronic res- 
onance gas(HRG) EOS comprising all the resonances be- 
low the mass 2.5 GeV. The entropy density of complete 
EOS is obtained as 



s = 0.5[1 + tanh(x)]sHRG + 0.5 [1 — tanh(x)]s lattice 

(6) 

with x= T ^ c , AT = 0.1T C . Compared to lattice simu- 
lation, entropy density in HRG drops slowly at low tem- 
perature, trace anomaly e ~4 P drops faster in lattice sim- 
ulation than in HRG model. It is difficult to resolve 
whether the discrepancy between lattice simulation at 
low temperature and HRG model is due to failure of HRG 
model at lower temperature or due to the difficulty in re- 
solving low energy hadron spectrum on a rather coarse 
lattice [|| . 




TABLE I: Initial central energy density (sj) and temperature 
(Ti) of the fluid in b=0 Au+Au collisions, for different values 
of viscosity to entropy ratio (r)/s). The predicted <f> meson 
multiplicity and mean pr are also noted. They should be 
compared with STAR measurements, ^j— = 7.95 ± 0.74 

ay ex 

and < p% > ex = 0.977 ± 0.064. 



rj/s 





0.08 


0.12 


0.16 


ei(GeV/fm 3 ) 


35.5 


29.1 


25.6 


20.8 




± 5.0 


± 3.6 


± 4.0 


± 2.7 


Ti (MeV) 


377.0 


359.1 


348.0 


330.5 




±13.7 


±11.5 


±14.3 


±11.3 


dN<f 
dy 


7.96 


8.01 


8.22 


8.13 


< pT<t> > 


1.019 


1.062 


1.111 


1.174 



C. Initial conditions 

Boost-invariant solution of EqsH] and [2] require initial 
conditions, e.g. transverse profile of the energy density 
(e{x, y)) and fluid four velocity (u(x, y)) and stress tensor 
ir^(x,y) at the initial time (t^). Relaxation equation 
(Eq|2]) require to specify the relaxation time . A freeze- 
out prescription, e.g. freeze-out temperature Tp is also 
needed. In the present paper, we fix the initial condition 
of the fluid as it was obtained in our analysis of <f> mesons 
At the initial time Ti=0.6 fm, the fluid velocity is 
zero, v x {x,y) — v y (x,y) — 0, the energy density of the 
fluid is distributed as, 

e(b, x, y) = £i[0.75N part (b, x, y) + 0.25N coU (b, x, y)], 

(7) 

where N par t(h, x, y) and N co a (b, x, y) are transverse pro- 
file of the participant and collision number distribution 
respectively, in an impact parameter b Au±Au collision, 
calculated in a Glauber model. Si is the central energy 
density in b — collisions. The shear stress tensor is 
initialised to boost-invariant value. For the relaxation 
time we use Boltzmann estimate t„ = Zrj/Ap. Freeze-out 
temperature is chosen to be Tp=150 MeV. The central 
energy density Si is obtained by fitting <f> multiplicity in 
0-5% Au±Au collisions The fitted values of central 
energy density and temperature are noted in tablelH As 
expected, the central energy density or temperature is 
reduced in more viscous fluid. 



III. RESULTS 

With the initial conditions as described above, we have 
computed transverse momentum spectra of pions, kaons, 
protons and 4> mesons. In Fig[3l predicted pion spectra 
from ideal fluid evolution in 0-5% Au±Au collisions are 
shown. Thermal pion's and decay pions are shown sepa- 
rately. Decay pions contribute mainly at low px < lGeV. 
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FIG. 5: (color online) same as in Fig|4]but for kaons. 



Decay and thermal pions together (the solid line) reason- 
bly well explain the PHENIX data [l(| ■ Note that initial 
condition of the fluid was not tuned for 7r mesons, Rather 
it was tuned to fit <p meson multiplicity in 0-5% collisions. 

In FigQ] present model predictions for ff "^ 7r , in 
0-60% Au+Au collisions are compared against the 
PHENIX data [l(| • For the ease of computation and with 
the understanding that resonances contribute mainly at 
low pt < 1 GeV, we have omitted decay pions. Accord- 
ingly spectra are shown only in the pt range 1 < pt < 3 
GeV. In FigH the colored symbols are the PHENIX 
data The solid, dashed, medium dashed and long 

dashed lines are hydrodynamics model predictions for the 
spectra in ideal fluid and in viscous fluid with viscosity 
to entropy ratio ry/s=0.08, 0.12 and 0.16 respectively. 
At large momentum pion yield increases with viscosity. 
From FiglU it is evident that in central and mid-central 
collisions, both the ideal and minimally viscous fluid rea- 
sonably well explains the data. Data are over predicted 
in fluid evolution with i]/s=0.l2 and 0.16. To obtain a 
quantative idea of fit to the data in ideal and viscous hy- 
drodynamics, we have computed x 2 /A?\They are noted 
in table [TTJ For all the collisions centrality, minimum 
X 2 /N is obtained in fluid evolution with the ADS/CFT 
limit of viscosity rj/s= 0.08. x 1 /N is comparatively large 
in ideal fluid. Compared to ADS/CFT limit of viscos- 
ity, fit to data is considerably poor evolution with vis- 
cosity, 7//,s=0.12 and 0.16. We also note that irrespec- 
tive of viscosity, quality of fit gets poorer as the colli- 
sions become more and more peripheral. For example, 
in minimally viscous fluid, upto collision centrality 30- 
40%, x 2 / N ~ 10, but beyond 30-40% collision centrality, 
X 2 /N increases by a factor of ~2. 30-40% Au+Au colli- 
sions corresponds to average impact parameter b=8.3 fm. 
It appear that in a hydrodynamic model, pion spectra is 



TABLE II: The x 2 /W values, in ideal and viscous (2=0.08- 
0.16) evolution for the it, K, p and cj> in different centrality 
ranges of Au+Au collisions are shown. In the last two rows, 
X 2 /N values for the combined data sets, (ir+K+p+cj)) and 
(tt+ K+ (f>) in the centrality range 0-40% are shown. 
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explained only in b <8.3 fm Au+Au collisions. 

Compared to pions, resonance production of kaon is 
much less. In FigEJ PHENIX data [Icj for K+ % K " are 
compared with hydrodynamic model predictions. Fit to 
the data in the entire pt range are shown. Kaon spectra 
in the centrality range 0-40% are also well explained in 
ideal and minimally viscous fluid evolution. As it was 
for pions, in more viscous fluid {rj/s >0.12) description 
deteriorates. Again, to obtain a quantative idea of fit 
to the data, we have computed x 2 /N. To be consistent, 
with the analysis of pion spectra, in the \ 2 analysis, data 
only in the pr range 1 < pt < 3 GeV are included. Val- 
ues are noted in tableHl Unlike for pions, best fit to the 
kaon data is obtained in ideal hydrodynamics. Viscous 
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FIG. 6: (color online) same as in Fig|4]but for protons. 



hydrodynamic gave comparatively poorer fit. 

In Fig|6l hydrodynamic model prediction for proton 
spectra are compared with the PHENIX data pi| for 
Interestingly, both ideal or viscous hydrodynamics, 
with the lattice based EOS fails to reproduce the proton 
data. In central or mid-central collisions, data are under- 
predicted by a factor of ~2. Proton spectra in peripheral 
collisions are comparatively better fitted. It is reflected in 
the X 2 /N analysis also (see table JTTJ) . Minimum x 2 /N for 
the proton spectra is order of magnitude larger than that 
obtained for pion or kaon spectra. We have neglected res- 
onance production. However, resonance contribution to 
proton is very small and even if resonances are included, 
model predictions will not agree with the experiment. In- 
ability of the model to correctly predict proton spectra is 
possibly due to the neglect of baryons in the model. We 
have assumed a baryon free fluid. However, the matter 
produced in Au+Au collisions at RHIC are not entirely 
baryon free. Note that with baryons in the fluid, proton 
production will depend exponentially on the chemical po- 
tential. Even a small chemical potential at the freeze-out 
will considerably enhance proton production. Inability of 
the model to correctly predict pr-spectra of proton may 
also be due to the equation of state used presently. The 
equation of state is a combination of a recent lattice simu- 
lation and non interacting hadron resonance gas. The 
lattice simulation was performed with two light quark 
and a heavier strange quark. Light quark masses that 
are nearly twice the physical masses, the strange quark 
mass was adjusted to its physical value. Consequently, 
pion is heavy, «220 MeV and kaon mass is physical 
M K «503 MeV 0. Though it is not mentioned in [2j, we 
expects that protons are also heavier than the physical 
ones. Then in the lattice based EOS the confined phase 
is a fluid of heavy pions and protons but physical kaons. 
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FIG. 7: (color online) same as in Fig(4]but for (f> mesons. 



At the freeze-out fluid cells will contain less number of 
heavier protons than it would have otherwise and proton 
production will be reduced. Pion production will also be 
reduced, however, being a lighter particle, the effect of 
heavy pion mass will be less pronounced. For example, 
if we approximate dN/dycPpx oc exp(— \J m 2 + p T /T) , 
then in the pr range 1-3 GeV, for ~ 50% increase in 
pion mass, production is reduced only by 3-10%. For a 
similar increase in proton mass production is reduced by 
60-90%. Indeed, even if proton is only 20% heavier than 
physical proton, proton production is reduced by 30-60% 
in the pr range 1-3 GeV. More detailed study is needed 
to understand the proton spectra. 

As mentioned earlier, we have initialised the fluid (ideal 
or viscous) to reproduce <f> meson multiplicity in 0-5% 
Au+Au collisions. In FigJTl cj) meson's px spectra are 
studied. As before, the solid, dashed, medium dashed 
and short dashed lines are from evolution of fluid with 
J7/s=0,0.08, 0.12 and 0.16 respectively. The filled cir- 
cles are from the STAR experiment [Tlj . For the <j> data 
also, pt spectra are best explained in ideal fluid evolu- 
tion. Comparatively poor description is obtained in vis- 
cous evolution. It is evident also from the x 2 values in 
tableHIl In all the collision centrality, x 2 /N is minimum 
in ideal fluid evolution. % 2 /7V's are comparatively larger 
in minimamly viscous fluid evolution. For ?y/s=0.12 or 
0.16, compared to ideal fluid, in viscous fluid evolution, 
X 2 /N increases by a factor of 3-6. 

In FiglHJ we have shown the x 2 /A values for the com- 
bined data sets; w, K , p and <f>. Collision centrality upto 
30-40% are included only. As indicated earlier, hydrody- 
namic description to the data gets poorer beyond this col- 
lision centrality. In Fig 13 the filled square are the % 2 /A 
of the combined data sets, as a function of viscosity rj/s. 
X 2 /N analysis definitely indicate that identified parti- 
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FIG. 8: Filled squares are the \ 2 /N for the combined data 
sets (-7T, K, <j> and p) as a function of viscosity to entropy 
ratio (t)/s). The filled circles are the same when proton data 
are excluded. Experimental data only upto 30-40% collision 
centrality are included. 



cle Pt spectra do not demand large viscosity. Minimum 
X 2 /N 32 is obtained in minimally viscous evolution. 
But comparable description to the data is also obtained 
in ideal fluid evolution. In FigEl the filled circles are the 
X 2 /N when proton data are excluded from the analysis. 
As noted earlier, proton spectra are not well reproduced 
in the model. If proton data are excluded, x 2 /N values 
improves. Again the best fit to the combined tt, K and 
(j> data is obtained in minimally viscous fluid evolution, 
X 2 /N ~ 19. Ideal hydrodynamics give comparable fit. 
The results are consistent with recent estimate of QGP 
viscosity [3]. In [ij analysing <f> meson data, it was con- 
cluded that nearly perfect fluid is produced in Au+Au 
collisions at RHIC energy. Transverse momentum spec- 
tra of identified particles also lead to similar conclusions, 
in Au+Au collisions, a nearly perfect fluid is produced. 

Before we summarise our results, it is important to 
mention that we have neglected bulk viscosity. Exper- 
imental data, which include the effect of bulk viscosity, 
if there is any. In general, bulk viscosity is an order of 
magnitude smaller than shear viscosity. But in QGP, it 
is possible that near the cross-over temperature, bulk vis- 
cosity is large [111, H3] ■ Effect of bulk viscosity on particle 
spectra and elliptic flow is studied in [25j . It appears that 
even if small, bulk viscosity can have visible effect on par- 
ticle spectra and elliptic flow. Neglect of bulk viscosity, 
will artificially increase the effect of (shear) viscosity. In 
other word, if bulk viscosity is included, comparable fits 
to the data can be obtained with still lower value of rj/s. 
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IV. SUMMARY AND CONCLUSIONS 

To summarise, in a hydrodynamical model, where the 
evolution is governed by a lattice based equation of state 
with a confinement-deconfinement cross-over at tempera- 
ture T co =196 MeV, we have analysed the transverse mo- 
mentum spectra of identified particles, e.g. pions, kaons, 
protons and 4> mesons. It is assumed that Au+Au col- 
lisions produce a 'baryon free' ideal/viscous fluid. Ideal 
or viscous (^/s=0.08-0.16) fluid was initialised to repro- 
duce <j) meson multiplicity in a central (0-5%) Au+Au 
collision. Hydrodynamic evolution of the ideal or min- 
imally viscous (77/5=0. 08) fluid, initialised to reproduce 
4> multiplicity in 0-5% Au+Au collisions, reasonably well 
reproduces transverse momentum spectra of tt, K and 
<f> in central and mid-central collisions. In peripheral 
collisions, 40-50% and beyond, the description to the 
data gets poorer. Description to the that data is also 
poor in evolution of fluid with viscosity larger than the 



ADS/CFT limit. Hydrodynamical evolution of baryon 
free ideal or viscous fluid however do not generate enough 
protons to agree with experiment. Proton spectra are 
underpredicted by a factor of 2. Poor fit to the pro- 
ton data is possibly due to the neglect of baryons in the 
model. Fluid produced in Au+Au collisions at y / s=200 
GeV is not entirely baryon free. It is expected that the 
fits to proton data will improve if baryons are included 
in the model. Poor fit to proton data may also be due to 
comparatively large light quark masses in lattice simula- 
tion. Light quarks are approximately twice the mass of 
physical quarks, consequently protons are heavy. More 
detailed study is need to sort out the issue. Our analy- 
sis also indicate that the transverse momentum spectra 
of the combined data set, (tt, K, p and 4>) or (tt, K, 
and (/>), in 0-40% collision centrality are best explained 
in minimally viscous fluid. Nearly equilvalent description 
is obtained in ideal fluid evolution. Data definitely reject 
large viscous fluid, -q/s >0.12. 
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